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Abstract 



We expect that s-channel unitarity should materiaUze in hard DIS reactions through 
screening corrections (SC) indicating that the gluon distribution function is approaching 
saturation, it is not as yet clear what the kinematical scales are at which these effects 
become important. While the global DIS 7*p total cross section, or F2{x,Q'^), data are 
well reproduced by DGLAP evolution without substantial SC, there exists experimental 
data from HERA which suggests deviations from DGLAP predictions in the small 
and X limits. These signatures are observed in both the fine details of F2{x, Q"^) provided 

and x are small enough, as well as in the diffractive channels. In this investigation 
we present a detailed study of 9^2/5 In which is supported by a coupled analysis of 
J/^* photo and DIS production. Both channels are directly proportional to xG{x,Q'^), 
and as such serve as excellent discriminators between different approaches and models. 
In the first phase of our investigation we have found that none of the latest editions of 
the parton distribution functions (GRV98, MRS99, CTEQ5) provides an adequate and 
simultaneous reproduction of logarithmic slope of F2 at small values as well as J/^' 
photoproduction (Details of this will be published separately We then show that 

taking GRV98NLO as input and correcting it for SC, we can reproduce the recent HERA 
data well. The calculation depends on one parameter = 8.5GeV~'^ which is directly 
determined from the J/^ photoproduction forward differential slope. With this input we 
obtain an excellent fit to the J/^ photo and DIS production data. Our calculations made 
in the LLA of pQCD take into account the corrections implied by the real part of the 
production amplitude, off diagonal (skewed) gluon distributions and the Fermi motion 
of the charm quarks within the bound Charmonium system. The SC are consistently 
calculated for both the percolation of a qq through the target and the screening of the 
gluon parton distribution which forms the base of our calculation. Our main conclusion 
is that, whereas we find strong support for the need for SC in the small and x limits 
of the channels we have investigated, the latest HERA data is not sufficiently precise to 
directly determine the gluon saturation scale. 



1 Introduction 



Over the past few years we have been witness to vigorous experimental, phenomenological and 
theoretical investigations of the proton deep inelastic scattering (DIS) structure functions and 
some exclusive channels with small < 5GeV^ and x < 10~^. These comprehensive studies 
aim at establishing the applicability and possible need for a re-formulation of pQCD, as we 
know it, when approaching the kinematic interface with the less understood npQCD dominated 
domain. The standard procedure for the pQCD analysis of DIS on a nucleon target has been 
to utilize the DGLAP evolution equations for the structure functions as the key ingredient for 
fixing the parton distribution functions (p.d.f.). These p.d.f. are then used as input for the 
calculations of exclusive DIS channels, usually executed in the color dipole approximation. 

The physics of small and small x is associated with the search for the scale of gluon 
saturation implied by s-channel unitarity [|]. One should remember, though, that gluon sat- 
uration signals the transition from perturbative to non perturbative QCD. This transition is 
preceded by SC signatures which are expected to be experimentally visible even though the 
relevant scattering amplitude has not yet reached the unitarity black limit. Moreover, from 
our experience with soft Pomeron physics, we know that different channels have different scales 
at which unitarity corrections become appreciable. Specifically, the scale associated with the 
diffractive channels are considerably smaller than those associated with the elastic channel 0] . 

Inspite of significant theoretical progress in recent years 0, it is still not clear what the 
saturation scale in the present experimentally accesible kinematic region is. We recall that, 
while the global analysis of F2{x,Q'^) (or (Jtof (W^ Q"^)) data shows no conclusive deviations 
from DGLAP, there are dedicated HERA investigations suggesting deviations from the DGLAP 
expectations in the small and x limits. These signatures are observed in both the fine details 
of F2{x, Q"^) as well as in the diffractive channels, provided and x are small enough. 

In the following we present a detailed study of dF2/d\nQ'^ which is supported by a coupled 
analysis of J/\l/ photo and DIS production. The strategy of our investigation is based on the 
observation that these observables ( in LLA of pQCD ) are linear in xG{x, Q'^) and {xG{x, Q^))^ 
respectively, and being relatively well measured may serve as effective discriminators when we 
compare their detailed features with existing relevant theoretical approaches and models. In 
the first phase of our investigation, we found that none of the latest p.d.f. provides an 

adequate simultaneous reproduction, at small pQCD scales, of the recent HERA data on the 
logarithmic slope of -F2 @]0 at small and x, as well as the abundant high energy data on 
J/m photo and DIS production lTl|] 0. We then proceed to show that when GRV98NLO is 
corrected for SC, as suggested in our previous publications |]13||M, it gives a good description 
of these data. The SC calculation depends on one parameter, R^= S.SGel^"^, which is directly 
deduced from the J/'^ photoproduction forward differential slope. We elaborate on a recent 
suggestion ||T^ that the gluon saturation scale may be determined by examining the behaviour of 
dF2{x, Q"^) /dlnQ"^ against and x at fixed W values. In our opinion, this suggestion which was 
discussed also in a recent presentation of the ZEUS new data [^], reflects the particular kinematic 
relationship between , W"^ and x and does not provide an unambiguous determination of the 
desired scale. 

Our approach can be tested in the high energy analysis of photo and DIS exclusive produc- 
tion of J /'^ . This is seemingly a straight forward procedure as the LLA pQCD calculation of 

this cross section is proportional to {xG{x, Q^)) , where x and are determined from mc, the 
c-quark mass, and W, the cm. energy. A realistic calculation depends on a few corrections to 
the original pQCD estimate: C\ due to the amplitude real part, due to the contribution of 
the off diagonal (skewed) gluon distribution and due to the Fermi motion deviations 



from the simple static non relativistic estimate of the J /'^ wave function |17|. All of these cor- 
rections contribute to the final (amplitude squared) estimate. We present a detailed analysis 
of this channel which is consistent with our approach and input choices made in the analysis 
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of dF2/d\nQ^. 

We conclude with some general comments and observations on the outstanding problems of 
gluon saturation and screening corrections. 



2 The small and x behavior of dF2/d\nQ'^ 

Checking for unitarity corrections in F2 studies is not simple. As is well known, a global DGLAP 
analysis of the data with the recent p.d.f. is adequate. A study fl^ comparing the 

screened and non screened DGLAP calculations of F2{x, Q^) showed only a small difference due 
to SC even in the small and x attained by present HERA measurements. Clearly, a unitarity 
study in the above kinematic limit requires a dedicated investigation, confined to small and 
X, which can magnify the presumed experimental signatures. We recall that in the small x limit 
of DGLAP we have 

^^i:^ = ^xG^^^^^(x,Q^). (1) 
olnQ'^ 97r 

Accordingly, a significant deviation of the data from Eq. (|I]), where xG^'^^^^ is obtained 
from the global F2 analysis, may serve as an experimental signature indicating the growing 



importance of unitarity corrections. This was first suggested by Caldwell |T^, showing a rather 
complicated plot of dF2/dln in which each point had different and x values. The Caldwell 
plot suggested a dramatic turn over of dF2/dlnQ'^ corresponding to of about 3GeV^ and 
a; < 5 ■ 10~^ in contrast to the behavior expected from GRV94 at sufficiently small 
and X. The problem with this presentation is that, as suggestive as it may seem, it does 
not discriminate between different dynamical interpretations [|T^] |TB[ ||2T| . It is actually 
compatible with an overall data generator as well as the latest p.d.f. which were re-adjusted 



to account for this observation. We conclude that, where as a reproduction of the Cadwell plot 
is required as a pre condition for a serious consideration of any suggested model, the plot on its 
own can not serve as an effective discriminator between models even in the extreme case when 
they are fundamentally different. 

A far better discrimination is obtained if we carefully study the small and x dependences 
of dF2/ dlnQ"^ at either fixed or fixed x values being free from the kinematic correlation 
between and x that plagued the Caldwell plot. Such preliminary HERA data have recently 
became available [pjP]. As we shall show, a pQCD analysis of these data is consistent with a 
SC interpretation jl^ . 

We follow the eikonal SC formalism presented in Ref . ||14j| , where screening is calculated in 



both the quark sector, to account for the percolation of a qq through the target, and the gluon 
sector, to account for the screening of xG{x, Q^). The factorizable result that we obtain is 

SC in the quark section are given by 

= ^i,GO°"'"(x,Q2)r((,2). (4) 



The calculation is significantly simplified if we assume a Gaussian parameterization for the two 
gluon non perturbative form factor, 

m = ^e-^v«^ (5) 
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SC in the gluon sector are given by 

Q2) = d,{x, Q2)a;G'^^^^^(x, Q^), (6) 

where 

TV^ Jx X' Jo J ^ ' 

Note that Kg[x' .Q'"^]}?) = Q'^; 6^) defined in Eq. (|). An obvious difficulty in the 

above calculation of xG^^ stems from the fact that the Q'^ integration spans not only the 
short (pQCD), but also the long (npQCD) distances. To overcome this difficulty we assume 
that 

-i2 ^ ..2\ _ V_ 



xG (x, Q2 < /i^) = ^xG (x, Q^) , (8) 



where /z^ ~ Qg. Our choice of the above interpolation is motivated by the gauge invariance 
requirement that xG oc when —>■ 0. 

The SC calculation described above can be applied to any given input p.d.f. where the only 
adjusted parameters are and /i^. As we shall see in the next section, = 8.5GeV~^ is 
determined directly from the forward slope of J/\l/ photoproduction. fi"^ is conveniently fixed at 
Qq, the lowest value accessible for the input p.d.f. we use. Once we have chosen our p.d.f., 
our SC calculation is essentially parameter free. We have checked that our output results are 
not sensitive to the fine tuning of these fixed parameters. 

Our results are presented in Figs. 1,2 and 3. Throughout this investigation we have used 
as input the MS version of GRV98NLO. Using the GRV98DIS version provides very similar 
results. Following are some comments relating to our calculations and results: 
1) As can be seen, in the limit of small and x there is a significant difference between the 
screened and non screened values of dF2/ dlnQ'^ . As expected the SC results are smaller and 
softer than the non screened input. 2) Visibly, our overall reproduction of the experimental 
data is very good, in particular when considering that our input is essentially parameter free. 
A proper calculation requires the knowledge of the unknown theoretical errors. 3) If we 
follow the standard procedure and replace the theoretical errors with the experimental ones, we 
obtain excellent 'x^ /ndf = 0.75 for 21 HI data points with the exception of 3 points which are 
visibly out of line. The ZEUS data has errors which are considerably smaller and as a result 
our /ndf is not as good, even though our reproduction of the ZEUS data is reasonable. 
4) The ZEUS Q"^ = 1.9GeV^ and HI = 3.0GeV^ data are somewhat softer than our 
predictions, which do not contain a soft non perturbative background. 



Golec-Biernat and Wuesthoff [|T^] have suggested studying the and x behaviour of 
dF2/ dlnQ"^ at ffxed W as a method to determine the gluon saturation scales from the an- 
ticipated turn over in these plots. Recent ZEUS low presentations of these plots show, 
indeed, the anticipated turn over structure in these figures, seemingly suggesting that gluon 
saturation is attained at ~ IGeV"^. In our opinion the proper variables to study F2 are 
X and Q"^. Trying to study the structure function by introducing other variables, such as W, 
may result in spurious effects which are predominantly kinematic. In the specific procedure 
suggested by Golec-Biernat and Wuesthoff, the combination of the kinematic relation between 
X, and W with the very general behaviour of xG{x, Q"^) is sufficient to produce a turn over. 
Its exact location depends on the details of the numerical input. Consequently, the suggested 
fixed W plots do not convey any new dynamical information even if such information is hidden 
in the analyzed data. Actually, it seems that any F2 (or xG) parameterization which factorizes 
the and x dependences, such as the Buchmueller- Haidt model ^4] , is capable of producing 
the fixed W turn over effects. 

To illustrate this point we consider the fixed W behaviour in two models which have very 
different dynamics: 
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1) Our own GLMN , which is a pure pQCD dipole model with SC. As such, our model 
relates indirectly to gluon saturation, even though it is constructed so as to include unitarity 
corrections below actual saturation. 



2) The DL two Pomeron parametrization [^, which is based on the Regge formalism, and 



consists of the coherent sum of contributions from a "hard" and a "soft" Pomeron, a normal 
Reggeon and an additional contribution from the charmed sector which is proportional to 
the "hard" Pomeron. Each of these fixed j-poles are multiplied by a fitted form factor. 
The parameters of the "model" were determined from the requirement of a "best fit" to the 
experimental -F2(x, Q^) data, and are not directly associated with the gluon distribution in the 
proton or " saturation" . 

The results obtained for the logarithmic slope of F2 at fixed W for both parametrizations 
are compared to the ZEUS experimental results in Fig. 4 and 6. In Figs.5 and 7 we also display 
the behaviour of dF2/ dlnQ'^ at fixed and fixed x. We note that both "models" provide 
a reasonable description of the ZEUS data including the observed turn over in the fixed W 
plots. As GRV98 is only applicable for Ql > IGeV'^, we have repeated the GLMN calculation 
with GRV94 which has Ql = OAGeV^ and again reproduce the ZEUS fixed W turn over effect. 
From the above it appears that any model which provides a reasonable description of F2{x, Q^) 
will exhibit a turn over in fixed W plots of the logarithmic derivative of -^2(0;, Q"^) , this occurs due 
to the relation between the kinematic variables x, and W, and is not related to "saturation". 

However, the examples (the Buchmueller-Haidt model and Donnachie-Landshoff approach) 
which have been used to demonstrate that the "saturation" is not a unique mechanism for a 
turn over at fixed W, have one common feature: the "soft" contributions are concentrated at 
a rather large typical scale > 2GeV^. This observation supports the idea that the so called 
soft Pomeron stems from rather short distances ||25|| . 

Therefore, an equally good description of the experimental data can be achieved either from 
gluon "saturation" ( shadowing corrections ) or by matching soft and hard processes if the soft 
ones occur at rather short distances < 0.5 GeV~^ . 

In conclusion, whereas we find strong support for the need for SC in the small and x 
limits for dF2/ dlnQ"^ we are unable, as yet, to determine the gluon saturation scale directly 
from the latest HERA data. The gluon saturation scale may be theoretically estimated from 



the contours produced at the boundary of discussed in our papers [l^ 



3 Photo and DIS production of J/^ 

The t=0 differential cross section of photo and DIS production can be calculated in the dipole 
LLA approximation p6[| [p^ p8[ | p9|| ||30|| . When using the static non relativistic approximation 
of the vector meson wave function [|26|||28|, the differential cross section has a very simple form, 
viz. 



Jt 



Q 



where in the non relativistic limit we have 



M2 + Q2 



(9) 



(10) 



and 



X 



452 



In the following we discuss in some detail the photo and DIS production of J/\E'. There is 
a lot of data available for this channel [|lOl |Jll| for the integrated cross section spanning a 
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10-4 10-3 lQ-2 IQ-l lQ-4 lQ-3 IQ-B IQ-l 
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Figure 2: x dependence of ZEUS logarithmic slope data at fixed compared with our calcula- 
tions. 
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Figure 3: dependence of HI and ZEUS logarithmic slope data at fixed x compared with our 
calculations. 
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Figure 4: ZEUS logarithmic slope data at fixed W compared with our SC calculation. 
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Figure 6: ZEUS logarithmic slope data at fixed W compared with DL two Pomeron model. 
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relatively wide energy range. From a theoretical point of view, its hardness (or separation) scale 
is comparable to the scales we have studied in our dF2/ dlnQ'^ analysis. To relate the integrated 
cross section to Eq. (^ we need to know B, the J/\E' forward differential cross section slope. 
For this we may use the experimental values, which are approximately constant. As we shall 
see, SC account well for the reported moderate energy dependence [|10| . 



The main problem with the theoretical analysis of J/'^ is the observation that the simple 
pQCD calculation needs to be corrected for the following reasons: 

1) A correction for the contribution of the real part of the production amplitude. This correction 
is well understood and is given by C\ = (1 + p^), where p = ReA/ImA = tg{^) and 
A = dln{xG)/dln{l). 



2) A correction for the contribution of the skewed (off diagonal) gluon distributions |[T6|. This 
correction is calculated to be 

- [ 0Fr(A + 4) ) ■ ^^2) 

3) A more controversial issue relates to the non relativistic approximation assumed for the J/^ 
Charmonium. Relativistic effects produced by the Fermi motion of the bound quarks result in 
a considerable reduction of the calculated pQCD cross section [|17| . We denote this correction 



Cp and note that it is very sensitive to the value of m^. Ref . IT? | assumes that ~ 1.5GeV 
and obtains Cp ~ 0.25 with minimal energy dependence. A small change in the input value 
of rric changes the above estimate significantly. We suggest, therefore, to consider Cp as a free 
parameter. In our calculations we have used Cp = 0.66 which corresponds to a c-quark mass 
of approximately 1.53GeV. 

Since the J/ /"^ photo and DIS cross sections are proportional to {xG{x,Q'^))'^ , the study 
of this channel can serve as a compatibility check supplementing our study of dF2/ dlnQ"^ . 
As noted in the introduction our SC approach was triggered by the observation that none of 
the latest p.d.f.'s can provide a good simultaneous reproduction of the two channels under 
consideration. 

Our calculation of SC for J/\l' photo and DIS production is rather similar to the loga- 
rithmic slope calculation presented in the previous section. We follow our earlier publication 



3T| and define the damping factors due to the screening in the quark sector i.e. the percolation 
of the cc through the target. This is given by the following expressions for the longitudinal and 
transverse dampings 



= ^ ' (13) 

and 



1 + (1 - ^)E,{^)e^. 



2 



^ T = ^ (14) 



Our expression for D^, the damping in the gluon sector is the square of the gluon damping 
defined in the previous section. 

Our final expression for the forward cross section is 

where Dq denotes the L and T components as appropriate. 
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Our calculations as compared with the data are presented in Figs. 8 and 9. As can be seen, 
our reproduction of the data is excellent with a 'x^ jn.d.f . which is well below 1. These excellent 

values are maintained when calculating over the entire data base as well as limiting ourselves 
to the high energy HERA data. Note that i?^ = S.SGel^"^, which is the essential parameter 
in the SC calculation is determined directly from the J photoproduction forward slope. In 
a model such as ours, we expect a weak dependence of Bu on the energy [^( see Fig. 10 ). 

4 Discussion 

We have shown that NLO GRV98 when screened both for the percolation of a qq pair through 
the target, and for multigluon exchange provides an excellent repoduction of both the HI 
and ZEUS data for the logarithmic derivative of F2{x,Q^) structure function of the proton. 
The experimental data and our model are consistent with dF2/ dlnQ"^ ( at fixed ) being a 
monotonic increasing function of 1/x. No deviation of this behaviour has been seen even at the 
lowest values of and x. 

The suggested turn over seen in dF2/dlnQ'^ at fixed W, as a function of x or Q^, does not 
appear to be connected with saturation effects, and cannot be used as a discriminator, as it 
appears in all parametrizations of F2 which provide a reasonable description of the data. 

The same screened NLO GRV98 also provides an excellent description of the photo and DIS 
production of J/^E', once corrections are made for the real part of the production amplitude, 
for the skewed ( off diagonal) gluon distribution and for the Fermi motion of the bound quarks 
in the charmonium wave function. 

In a seperate publication [jl| we show that we are unable to obtain a simultaneous good 
description of dF2/ dlnQ'^ and the photo and DIS production of J/'^ using other pdf's on the 
market e.g. MRS99, CTEQ4 and CTEQ5 even after including screening corrections. 

5 Conclusions 

Only screened NLO GRV98 is able to provide a satisfactory simultaneous description of the 
latest HERA data available for dF2/ dlnQ^ and photo and DIS production of J/'^ . 

The moral of the paper is that only by a simultaneous analysis of all data sensitive to 
shadowing ( saturation ) effects, can one obtain a reliable estimate of their sizes. 
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